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i£ Markov Chain Monte Carlo is repeatedly used to analyze the properties of 

q-i intractable distributions in a convenient way. In this paper we derive condi- 

tions for geometric ergodicity of a general class of nonparametric stochastic 
t— i ■ volatility models with skewness driven by hidden Markov Chain with switch- 

■ ing. 
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1. Introduction 

X ' 

The asymmetry in the intertemporal relation between volatility and the 
stock return is the purpose of active study. In general, asymmetric effects in 
volatility denote, that the effects of positive returns on volatility are different 
from those of negative returns of a similar magnitude. Asymmetry is also 
sometimes referred to us as a negative relation between the ex-post volatility 
in the rate of returns on equity and the current value of the equity, i.e., the 
commonly investigated leverage effect. It is therefore necessary to stress, 
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that although leverage effect exhibits asymmetry, not all asymmetric effects 
display leverage. 

Two different approaches are of considerable interest, namely the condi- 
tional and stochastic volatility models. In the class of conditional volatility, 
ARCH specifications that have been developed to capture asymmetric ef- 
fects are: the Exponential GARCH (EGARCH) model of Nelson (1991) [l8| 
and the GJR model of Glosten, Jagannathan and Runkle (1993) [8j. In the 
manner depicted earlier, EGARCH model can describe leverage whereas the 
GJR model can capture asymmetric effects but not leverage. On the other 
hand, stochastic volatility (SV) models have become a natural alternative 
to time-varying volatility of the ARCH family. The asymmetric property 
within the stochastic volatility framework is established on a straightforward 
correlation between the innovations in both returns and volatility. Unlike 
ARCH type models, the SV models allow volatility dynamics to follow some 
latent stochastic processes based on volatility's unobservable nature. Ba- 
sic stochastic volatility process (BSV) was originally developed by Hull and 
White (1987) [l3| and Jacquier et al. (2004) [lij in continuous time framework 
as a mode to capture negative correlation between the innovation terms in 
the Black-Scholes option pricing formula. This basic version of the SV model 
assumes stationary error disturbances and uncorrelated Gaussian white noise 
processes (Harvey and Ruiz, 1993 10])- In empirical research, extensions of 
a simple discrete time model due to Taylor (1986)|23|| have been analyzed by 
Wiggins (1987) (26( and Harvey and Shephard (1996) [20] in order to accom- 
modate the direct correlation. This generalization, based on the immediate 
correlation between the innovation is known as the SV with leverage (SV-L) 
model (Asai and McAleer, 2006 Other asymmetric models were sug- 

gested by Danielsson (1994) (if, which was similar in the spirit to that of 
the EGARCH model. Nelson (1991) [181 used the absolute value function 
to capture the sign and magnitude of the previous value of normalized re- 
turns in accommodating asymmetric behaviour into an ARCH-type model. 
Danielsson (1994)0 used the absolute value function as in Nelson (1991) 



181 ]. but incorporated the observed return into the SV specification as it is 
not computationally straightforward in the SV framework to incorporate the 
normalized disturbances, in a SV with leverage and size effects (SV-LSE) 
model. So, Li and Lam (2002) [il[ considered a different type of threshold 
effects model in which the breaks in the constant and autoregressive param- 
eter in the SV equation depend on the signs of the previous returns. An 
alternative form of asymmetry can be based on threshold effects, as pro- 
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posed in Glosten, Jagannathan and Runkle (1992) [8] in the context of condi- 
tional volatility models. A variety of symmetric and asymmetric, univariate 
and multivariate, conditional and stochastic volatility models is analyzed in 
McAleer (2005) [16]. Excess kurtosis and skewness in the residuals of the 
basic SV model has triggered the use of either the mixtures of normals or 
the central t-distribution (Geweke, 1994 Q; Shephard and Pitt, 1997 [l9(] ) - 
Hansen (1994) has considered skewness in a GARCH model using skew 
t-distribution errors allowing for both skewness and heavy tails to co-exist in 
a time- varying volatility specification. Also, Cappuccio et al. (2004) |4J have 
managed to generate skewness and excess kurtosis in the conditional distribu- 
tion of returns by assuming a skew-GED distribution. Volatility asymmetry, 
initially interpreted as the leverage effect, is expressed by the correlation of 
the returns with the future volatility. Under this assumption, when mean 
returns are negatively correlated with contemporaneous volatility, the nega- 
tive (positive) return shocks are associated with increased (decreased) future 
volatility. Harvey and Shephard (1996) [io[ have developed the asymmetric 
SV (ASV-t) model, which guarantees the existence of the martingale prop- 
erty. Recently, Wang et al. (2011) [25J have developed a heavy-tailed SV 
model with leverage effect, where a bivariate t-distribution, in the form of 
a scale mixture of normal, is used to model the error innovations of the 
return and volatility equations. Tsiotas (2012) j^] introduced generalised 
ASV models that incorporate financial datas stylized facts including excess 
kurtosis, skewness, and volatility asymmetry, using three alternative distri- 
butions: the noncentral-t (NCT), the skew-normal (SN), and the skew-t (ST) 
distributions. These ASV models incorporate both the Gaussian and the non- 
Gaussian assumption in the same specification, while taking into account the 
leverage effect. The asymmetric central-t distributed SV (hereafter ASV-t) 
model nests in the asymmetric noncentral-t distributed SV (hereafter ASV- 
nct) model. The other one is the asymmetric skewnormal SV (hereafter 
ASV-sn) model. Finally, in the asymmetric skew-t SV (hereafter ASV-st) 
model both excess kurtosis and skewness co-exist. The first multivariate SV 
model was proposed by Harvey, Ruiz and Shephard (1994) [llj, who specified 
the model in terms of instantaneous correlations in the mean and volatility 
equations. However, their estimation technique was based on the inefficient 
quasi-maximum likelihood (QML) procedure. Danielsson (1998) jgj suggested 
a multivariate SV-L model based on the specification considered by Harvey, 
Ruiz and Shephard (1994 nil), but only estimated a symmetric version of the 
model. Shephard (1996) j20| proposed a one factor multivariate SV model. 
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Parametric and log-volatility estimation of the asymmetric stochastic 
volatility models is usually implemented in a Markov Chain Monte Carlo 
(MCMC) set-up. By overcoming the problem of likelihood non-existence in a 
close form solution, because of normality departures in the observed process, 
the Metropolis-Hastings algorithm is common in the MCMC egine. These al- 
gorithms provide a reliable structure to explore the intractable distributions. 
In any MCMC analysis, the convergence rate of the associated Markov Chain 
is of practical and theoretical importance. Geometric ergodicity is essential 
in convergence of Markov Chains to their stationary distributions. A geomet- 
rically ergodic chain converges to its target distribution at a geometric rate. 
In addition to ensuring the rapid convergence required for useful simulation, 
geometric ergodicity is a key sufficient condition for the existence of Central 
Limit Theorems and consistent estimators of Monte Carlo standard errors, 

1. e., geometric ergodicity justifies the applicability of the Central Limit The- 
orems to ergodic averages along the path of the chain. The main aim of this 
paper is to provide a formal proof of geometric ergodicity for general non- 
parametric model, that incorporates both the asymmetry and is the natural 
extension to CHARME models provided by Stockis et al. 2010. 

2. The Skewed Stochastic Volatility-CHARME model 

Consider first a nonlinear pth order autoregressive process defined by 

X t = m{X t - U X t _ p ) + a(X t _i, X t _ p )e t + A{X t - U X t ^ p )i 2 t , (1) 

where m, a and A are unknown, real-valued Borel measurable functions. The 
e t and t t are mutually independent identically distributed zero-mean innova- 
tions. Since, it is unremarkably not realistic to postulate that the observed 
process has the same trend function m, volatility function a or skewness 
function A at each time instant, we consider a class of nonparametric time- 
series models, where between random change points, the process we observe 
is piecewise stationary. The dynamics of {X t } is driven by a hidden Markov 
chain {Qt} with values in a finite set {1,2, ...,K}. Our model is defined as 
follows: 

K 

-A* = ^ S tk (m k (X t _i, X t _ p ) + a k (X t _i, X t _ p )e t + A k (X t _ 1: X t _ p )tf) , 

(2) 
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where 

s _ f l,forQ t = k 
\ 0, otherwise 

m>k, o~k and = 1,2, ...,K are unknown functions, e 4 and i t are 

independent and identically distributed random variables with mean zero 
and they are mutually independent. 

2.1. Model definition 

We make the following assumptions, adopting Stockis et al. 22J point of 
view: 

Assumption 2.1. The process {Qt} is a first - order strictly stationary 
Markov chain which is irreducible and aperiodic with probability distribution 
(jTx, 7r 2 , ttk) and transition probability matrix A = {aij}i<ij<x- 

Assumption 2.2. Let Gt-i — o-{X s ,s < t — 1} be the a-algebra generated 
by {X s ,s < t — 1} and let G t -i G &t-i- Then P(Q t = j\Qt-i = i, Gt-i) — 
P(Qt = j\Qt-i — i) f or a M hj- I n other words, the hidden process Q t is 
independent of the past observations of {X t } given its own past. 

Assumption 2.3. Qt is uncorrelated with the et and it, given (Qt-i, X t -\, X t -2, 

Assumption 2.4. Both e t and i t are independent of X t _i, X t _ 2 , 

Assumption 2.5. The functions mj., and are bounded on compact 
sets for all k. 

Assumption 2.6. The iid random variables et have a continuous and every- 
where positive density f and the iid random variables i t have a continuous 
and everywhere positive density g. 

Without the loss of generality, we restrict ourselves to the case p = 1, i.e., 
rrik, o~k and Aj, are functions on the real line. We assume that Aj,(x) > 0. 
We should make the following assumptions 

Assumption 2.7. The iid random variables et and it have mean zero and 
variance equal to 1. Moreover, there exists the fourth moment E(if) = k. 
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Assumption 2.8. 



max lllil blip 1 

ie{i,2,...,i<r} \x\-++oc x 2 

Let <S( = (Sn, S t2 , S t K) T ■ We conclude that, under assumptions 1-4, the 
process 

Zt = X t ) T , 

which represents the transformed mixture process, is a Markov chain as well. 



Theorem 2.1. Under Assumptions 1-8, the process {Z t } is geometrically 
ergodic. 

Proof. We shall prove that the conditions of Theorem 15.0.1 (iii) of Meyn 



and Tweedie [17| are satisfied. 

1. The process Z t is (^-irreducible if we take ip as a product of the stationary 
probability distribution measure of {Qt} on {1,2, K} and the Lebesgue 
measure on R. Let A = Ai x A 2 be such that (p(A) > 0. A\ contains at least 
one integer k between 1 and K. It suffices to prove that there exists t such 
that for all k and 1 

P((S t+u X t+l ) e {e k } x A 2 \S! = e h X l = x) > 

with e k denoting a unit vector with the kth component equal to 1. 
We obtain 

P((S 2 ,X 2 ) e {e k } x A 2 \S X = e l ,X 1 = x) = 
P {Q 2 = k 1 X 2 eA 2 \S 1 = e ll X 1 = x) = 
P (X 2 eA 2 \Q 2 = k,S 1 = e l ,X 1 = x)P (Q 2 = k\Q 1 = I, X 1 = x) = 
aikP {m k (x) + a k (x)e 2 + A k (x)L 2 2 G A 2 ) = 



1 ( u-y-m k {x) \ 9{^A k (x))+9{ yj A&) 
aik / / — 7-r/ 7-v f , v I {y>0 }dydu 



A k (x) 

aikh k (x), 
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where hk(x) > 0. 

The one but last equality we obtain, because the density of m,k(x)+ak(x)e2 

is 

f(y) = —f ( y ~ mk{xY 

and the density of Ak(x)t% is 

= 7 / \ J {y>0}- 



Then we use a well-known formula for convolution of the distributions, i.e., 

(7 * g) (y) = f(u- y)g(y)du. 
Jr 

Similarly, 

P((S 3 ,X 3 ) e {e k } x A 2 \S X = e h X 1 = x) = 

P(Q 3 = k,X 3 eA 2 \S 1 = e l ,X 1 = x) = 
P (X 3 eA 2 \Q 3 = k,Q 1 = I, X 1 = x)P (Q 3 = fc|Qi = I, X, = x) = 

^2 ai i a i k \ \ p ( x > d y\® 3 = k iQi= j)P(y, du\Qs = k,Qi = j) = 

K 

aijajkhjkjx) 

3=1 

where 




and hj k (x) > 0. 
Finally, we obtain 

K 

P {{S t +i, X t +i) e {e k } x A 2 \Si = e h Xi = x) = aih—aj t -ikhj 1 ,...j t _ 1 (x) 

h,—k-i 

which is strictly positive for some t because of the irreducibility of {Qt} and 
the fact that ^ii,...j t _ 1 (^) > 0. 

2. Analogously can be proven the aperiodicity of {Z t }. 

3. Irreducibility and aperiodicity ensure that every petite set is small. It 
can be shown that in the model we consider, each compact set is small. It 
suffices to prove that for every compact set B such that (p(B) > and for 
every bounded Borel set A = A% x A 2 with <p(A) > there exists t such that 

ini [P((S t+1 ,X t+1 ) e {e k } x A 2 \S! = e l ,X 1 = x) > 0. 
Proceeding as in the proof of irreducibility, we need to show that 

K 

xeB * — ' 

iivii-i 

for some t. Continuity of hj lt j tl (x)a.nd irreducibility of {Qt} ensure us 
that the last infimum is positive (see @, S])- 

4. We will apply drift criterion of Theorem 15.0.1 (iii). We need to show that 
there exist L > 0, (3 > and a function V(Z) > 1 such that for ||Zt-i|| > L 
we have 

E{V{Z t )\Z t _ l = {e u x))-V{e u x) < 
V(e h x) 

Let V(Z t ) = l + X?.We obtain 

E{V{Z t )\Z t _ 1 = {e l ,x))-V{e l ,x) 
V{e h x) 

Ef=i ( m l(x) + + kA\{x) + 2m k {x)A k {x)) E (S tk \ 5 t _i = e z ) - x 2 



< 



1 + x 2 

Ef=i (mj(x) + a 2 k (x) + K^(ar) + 2m fc (x)A fc (3:)) a lk _ j 
The conclusion is obtained by Assumption 8. □ 



8 



3. Summary and remarks 



We have derived sufficient conditions for geometric ergodicity of a gen- 
eral class of asymmetric nonparametric stochastic processes with stochastic 
volatility. It is natural to ask, whether other kinds of ergodicity (eg. polyno- 
mial ergodicity) can also be related to the skewed SV-CHARME model. It 
is the purpose of our future work. 
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